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Abstract 

Let F{m\ rii , n2) denote the number of lattice walks from (0, 0) to (ni ,712), always staying in 
the first quadrant {(ni, 712); ni > 0, 712 > 0} and having exactly m steps, each of which belongs 
^ ■ to the set {E = (1,0),M^ = (-1,0), NE= {1,1), SW = (-1,-1)}. Ira Gessel conjectured that 

F{2n; 0,0) = 16" ^.^J?^ /J" . We pose similar conjectures for some other values of (711,712), 



o 
u 



(2)„(5/3)„ 

and give closed-form formulas for F(ni; m, 712) when 7ii > 712 as well as for F(2n2 — ni; 77i, 772) 
when 77i < 772- In the main part of the paper, we derive a functional equation satisfied by the 
generating function of F{m; 771,772), use the kernel method to turn it into an infinite lower- 
,S^ • triangular system of linear equations satisfied by the values of F{m; 77i,0) and F{m; 0,772) + 

j^ I F(rn; 0, 772 — 1), and express these values explicitly as determinants of lower-Hessenberg matrices 

with unit superdiagonals whose non-zero entries are products of two binomial coefficients. 

1 Introduction 
> 

^\ ' Let F{m; 771,772) denote the number of lattice walks from (0,0) to (711,712), always staying in the 

CN ■ first quadrant {(ni,n2); ni > 0,n2 > 0} and having exactly m steps, each of which belongs to the 

^ ; set {E = {1,0),W = (-1,0), NE = (1,1), SW = (-1,-1)}. From the obvious recurrence 

O, F(m; 711,712) = F(7?T, - 1; Til -I- 1,712) + F(7n - 1; 711 - 1,712) 

O' + F(m-1; 711 -Fl, 712 -M) -h F(7ii-1; 711 -1,712-1) (1) 






valid for rii > 1, 711,712 > 0, the initial conditions 



C^ ■ jprn . J 1, Til = 712 = 0, , . 

^(°'"i'"^) = jo, otherwise ^'^ 

(so called because we'll think of m as being the time variable), and the boundary conditions 

F(77i; 711,712) = 0, for tii < or 112 < 0, (3) 

we can calculate many values of F(r7i; 111,712). For example, the sequence F{2n; 0, 0)^o of the 
numbers of lattice walks returning to (0, 0) after 2n steps starts out as 

1, 2, 11, 85, 782, 8004, 88044, 1020162, 12294260, 152787976, 1946310467, 25302036071, 
334560525538, 4488007049900, 60955295750460, 836838395382645, 11597595644244186, ... 



Based on such empirical evidence, Ira Gessel conjectures (cf. p^, p. 4]) that 
Therefore we wih cah the numbers F{2n; 0, 0) Gessel numbers. 

2 Similar conjectures 

Analogous conjectures can be made about other points. For example, we conjecture that 
j,(.„. n n - i«n(l/2)„ / 5 (7/6), (llln2 + 183n-50) (5/6) \ 

More generally, looking at the points (0, k) it seems that 

F(2n;0,k) = 16" i^^^^" ( ,, /^^^V, ^ Pkin)+ ,, /^^^V, s Qkin)] 
^ ^ {k + 2)^\ii3k + 4)/3)/''^ ' ((3fc + 5)/3)/'^ ^y 

where Pk{n) is a polynomial of degree 2k — 2 and qk{n) is a polynomial of degree 2k. These two 
polynomial sequences seem to be non-holonomic. 
At (2n + 2k; 0, n) we seem to have 

i'{zn,u,n) - 4 (2n+i)(2)„ - (2)„ ' 

F(2n + 2;0,n) = ^^""H^+i)(3/2)„. ^ 

F(2n + 4;0,n) = ^"("+i)(s>|32n+33)(3/2)„ ^ 

77('o c. n ^ — 4"-i(n+l)(64n4+672n3+2648n2+4641n+3060)(3/2)„ 

r [2n -\- b; <j,ti) — WTE) ' ^^'^■i 

from which we conjecture that 

F{2n + 2k;0,n) = 4" .f(^^" rk{n) (5) 

where rfc(n) is a polynomial of degree 2k — 1 divisible by n + 1 for k > 1, and ro(n) = l/(2n + 1). 

It seems that this polynomial sequence is non-holonomic. 

Another empirical observation is that g{n) = F{2n + 1; 1,0) likely satisfies the second-order 

recurrence 

(n + 3)(3n + 7)(3n + 8) g{n + 1) 

- 8(2n + 3)(18n2 + 54n + 35) c/(n) (6) 

+ 256n(3n + l)(3n + 2) c/(n-l) = 0. 

According to algorithm Hyper, this recurrence has no hyper geometric solutions. But F{2n; 2, 0) 
seems to be non-holonomic. 

On the other hand, we seem to have 

F(n;n,0) = 1, 

F(n + 2;n,0) = i(n + l)(n + 4), 

F(n + 4; n, 0) = ^ (^ + 1) (^^ + ISn^ + 74n + 132) , 

F{n + 6; n, 0) = ^{n + 1) (n^ + 32n^ + 407n3 + 2620n2 + 8604n + 12240) , etc.. 



from which we conjecture that 

F{n + 2k; n,0) = Sk{n) (7) 

where Sk{n) is a polynomial of degree 2k with leading coefficient ^ttx+itT' which is divisible by n + 1 
when A: > 1. This polynomial sequence seems to be non-holonomic. While Sfc(O) = F{2k; 0,0) is 
hypergeometric by Gessel's conjecture, and Sfc(l) seems to be holonomic of order two as per ([6]), 
the sequences Sfc(2), Sfc(3), ... all seem to be non-holonomic. On the other hand, the coefficient 
sequences [n^'']sfc(n), [n'^^~^Sk{n), [n^'^~^]sfc(n), ... all seem to be hypergeometric, again harboring 
a polynomial sequence of increasing degrees. 

3 Some values of the numbers F{m; 711,712) 

Proposition 1 F{m; ni,n2) 7^ only if 
(i) ?n, = ni (mod 2), 
(ii) ni < m, 
(iii) 71-2 < \{ni+m). 

Proof: As F{m; ni,n2) 7^ 0, there exists a walk w from (0,0) to (^1,71,2) having m steps. 
Assume that out of these m steps, a, b, c resp. d are E, W, NE resp. SW steps. Then 



a + b + c + d = m, 








a — h + c — d = ni. 








C — d = 71,2, 








(mod 2). Also, ni = a - 


-b + c- 


- d<a+b+c+d= 


■- m, and 

D 



n2 = c — d<a + c= {m + ni)/2. 

Theorem 1 Let s(ni,n2) be the length of a shortest walk w from (0,0) to (ni,n2). 
(i) // Til > 7T,2 then w uses E, NE steps only, and 

F(ni; 7Zi,n2) = 5(721,71,2) 

(ii) // 7T-1 < 712 then w uses W , NE steps only, and 

T?i'^ \ ( \ ni + l (2n2 - ni + 1\ 
F[2n2-ni;ni,n2) = s{ni,n2) = — 

2n2 -7Zi + l\^ 712 + 1 / 

Proof: Denote by a, b, c resp. d the numbers of E, W, NE resp. SW steps in w, and by m the 
length of w. 

(i) 71-1 ^ ^2- By Proposition [TJl^ii) , m > ni. The walk consisting of 722 NE steps followed by 
ni — n2 E steps ends at (711,112) and has length ni, so it is the shortest such walk. Hence 
5(111,112) = E{ni] 711,112). From a + b + c + d = m = ui = a — b + c — d \t follows that 
6 = d = 0, so ID uses E and NE steps only. From a + c = m = ni and c = 112 it follows that 
a = Hi — 712- Thus w must contain ni — n2 E steps and 712 NE steps. There is no restriction 
on the order of these steps, so there are ("^) such walks. 



(ii) 71-1 <n2: By Proposition [T]||iii), m > 2n2 — ?^l• The walk consisting of n2 NE steps followed 
by n2 — ni W steps ends at (rii, 71,2) and has length 2^2 — ni, so it is the shortest such walk. 
Hence s(ni,n2) = F{2n2 — ni; 111,71,2). In the inequality m = a + b + c + d> —a + b-\-c — d = 
2(c — d) — {a — h + c — d) = 2n2 — ni equality holds when a = d = 0, so w uses W and NE 
steps only. Conversely, any walk w using W and NE steps only satisfies —h + c = ni and 
c = 7i2, so 6 = 722 — Til and ra = b + c = 2n2 — ni, implying that u; is a shortest such walk. 
Thus to compute the number 5(721,712) of these walks it suffices to enumerate lattice walks 
from (0,0) to (711,712) staying in the octant < ?ii < 71,2 and using W and NE steps only. 
We have the recurrence 

s(ni,?i2) = s(ni + 1,112) + s(7ii - l,?i2 - 1) for 112 > ni + 1 > 1 
with boundary conditions 



s(-l,n2) 
5(77-1,111) 



0, for 112 > 0, 

1, for 7ii > 0. 



To this end, we transform the octant < iii < 112 into the first quadrant by the linear map 

/i: (1,1)^(1,0), (0,1)^(0,1). 
The matrix corresponding to 

(1,0)^(1,1), (0,1) 



M 



(0,1) 



in the standard basis of . 



d2 • 



IS 



M" 



so the matrix corresponding to /U is 



M 



1 





1 


1 


1 





1 


1 



Writing n = (111,112), define 

n(iii,ii2) = s(M~^n) = s(7ii,iii + 112). 

Then 5(711,112) = u{Mn) = 11(111,112 — tii), and u satisfies the recurrence 

ii(7ii,ii2) = 5(111,111+112) 

= 5(111 + 1, Hi + 112) + 5(711 - 1, 111 + 712 - 1) 

= li(71i + 1,112 - 1) +li(71i - 1,112), 

for ni > A 112 > 1, with boundary conditions 



ll(-l,ll2) 
li(71i,0) 



5(-l,n2 - 1) = 
5(111,711) = 1, 



0, for 712 ^ 1; 
for 7ii > 0. 



(8) 



(9) 
(10) 



Let 

f{x,y) = Yl u{ni,n2)x'''y''^ 

ni,n2>0 

be the generating function of u. From ([8]) - (llOp we obtain in the usual way the functional 
equation 

{x -x'^ -y)f{x,y) = x-yf{0,y) (11) 

which can be solved by the kernel method. Since 

X -x + y = \x I \x I , 

substituting x = x{y) = ~ ^^ ^ ^^ ([TT]) yields 

/(„,,) = £M . ^-f^ = c(,), 

y 2y 

the generating function of Catalan numbers. Hence 

x-yC{y) 1 C{y) 



f{x,y) 



x-x'^-y rj. i+Vi-% l-xC{y) 



2 



n\„,n2 



Following [?, p. 154], this can be expanded into 

/(.,,)= f;.".cter.«= Y. "_' + ' p'" + '" + ').■.., 

n ^o 2?T'2 + ni + 1 \ no / 

ni=0 ni,n2>0 z i i i y z / 

so we read off 

, , ni + 1 (2n2 + ni + l\ 

u{ni,n2) = ■ — 

2n2 + ni + 1 y 722 J 

and, finally, 

ni + 1 / 2n2 — ni + 1 



F{2n2 — m; ni,n2) = s(ni,n2) = n(ni,n2— ni) 



2n2 — ni + 1 \ 71.2 — ^1 / 



Corollary 1 For all n >0, we have 
(i) F{n-n,0) = l, 
(ii) F{2n; 0,n) = Cn, the n-th Catalan number. 

Proof: By Theorem[Hi), F{n; n, 0) = (T) = 1. By Theorem [T]^ii) , 

1 /2n + l\ 1 f2n\ 

F{2n;0,n) = -—-[ = = C„. 

2n + l\n+l/ n+ l\ n 



D 



4 A functional equation for the generating function 



Let 



y"2 



(12) 



G{x,y,z) = J2 F{m;ni,n2)x"'y"'^: 

m,ni,n2>0 

be the generating function of the numbers F{m; ni, 71,2). From ([I]) - ([3j) we obtain in the usual way 
the following functional equation satisfied by G{x,y,z): 

K{x,y,z)G{x,y,z) = x{l + z)G{x, 0, z) + xG{x,y, 0) - xG{x, 0,0) - y z. (13) 

Here the polynomial 



K{x,y,z) = x{l + z){l + y z) -yz 



(14) 



is called the kernel of equation p3j) . 

In order to simplify (jlSh . we introduce another generating function 



H{x,y,z) = K{x,y,z)G{x,y,z) +yz. 



Since 



H{x,0,z) = x{l + z)G{x,0,z), 
H{x,y,0) = xG{x,y,0), 
H{x,0,0) = xG{x,0,0), 

equation (fT3|) becomes 

H{x,y,z) = H{x,0,z) + H{x,y, 0) - H{x, 0,0). 

This is the functional equation that we will work with in the sequel. Write 

H{x,y,z) = J2 ^(m; ni,n2)x™y"iz"^ 

m,ni,n2>0 

It follows from ([IT]), ^ and ^ that 

F{m;ni,n2) = if nin2 / 0, 

F{m;ni,0) = F{m — 1; ni,0), 

F{m;0,n2) = F(m - 1; 0,^2) + F(m - 1; 0,n2 - 1). 

Note that F(2n; 0,0) = F(2n + 1; 0,0), and Gessel's conjecture (jl]) can be stated as 



(15) 
(16) 



(17) 



G(x,0,0) 



H{x,0,0) 



3^2 



5/6, 1/2, 1 
2,5/3 



iFi 



Wx^ 



-1/2, -1/6 
2/3 



16x^ 



2x' 



In analogy to ([5]) we conjecture that 



F{2n + 2k + l; 0,n) = i'"' 



:i/2), 



rk{n) 



{k + 2) 
where rk{n) is a polynomial of degree 2/c+l, and the sequence of polynomials ^^(n) is not holonomic. 



6 



5 The kernel method 



Equations ()13p resp. ()17p cannot be solved right away because they seem to contain other unknown 
functions beside the fuh generating functions G{x, y, z) resp. H{x, y, z) (the additional unknowns 
being just sections of G{x,y,z) resp. H{x,y,z), of course). To obtain more information, we look 
for roots of the kernel w.r.t. one of the variables which are power series in the remaining variables. 
Substituting such roots into (fT3]l resp. (fTTj) yields additional equations which are free of the term 
containing the full generating function G{x,y,z) resp. H{x,y,z). 

In our case, the kernel (fill) is linear in x, and quadratic in y and z. The roots of K{x, y,z) = 
w.r.t. y resp. z are not power series in x, z resp. x, y. But solving K{x, y, z) = for x yields 



yz ^ y f z yz 

{l + z){l + y'^z) l-y\l + z l + yz 



( \ y^ y '^ y^ \ /iq\ 



n=lfc=l 

which is a power series in y, z satisfying x{0, 0) = 0. So we can substitute it into (jlTh to obtain 

H{x{y,z),0,z) + H{x{y,z),y,0)-H{x{y,z),0,0) =yz (19) 

where the rational function x{y,z) is given in (jlSp . This does not help us find H{x,y,z) (or a 
non-trivial section of it) directly. However, (J19p does determine all the coefficients of -ff (x, y, 0) and 
-ff (x, 0, z), and hence also of H{x, y, z) and G(x, y, z). How could we exploit this? 

6 Gessel numbers as determinants 



Expand the left-hand side of (llOp into power series in y and z, and equate the coefficient of y'^z" 
to (except for the coefficient oiy z which is equated to 1). This yields the following infinite system 
of linear equations for the values of F{m; ni,n2) on the planes ni = and n2 = 0: 



-m \ / —m 



?Ti — Ti (mod 2) 

u-m-n, ( _ „+^-„, F{m; m, 0) (20) 

2 / V^ 2 / 

7Ti+7i]^=ii (mod 2) 

= < „' ,, . ' foralln, u>0. 

I U, otherwise. 

Note that both sums are finite since a binomial coefficient vanishes when its lower symbol is negative. 
So m < u and n2 < v — u/2 in the first sum, and m + rii < u in the second. To put this system 
into a more compact form, pack its unknowns into an infinite matrix [/(i, j)]f^-^Q defined by 

f{i,j) = I ^;^;0'^'-^|' ^^^■' (21) 

^ ^ 1 F{j;i-j,0), t>j, ^ ^ 



or graphically 
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Inverting this transformation, we clearly have 

F(m; 0, n2) = /(m,m + 722), 
F{m; ni,0) = f{m + ni,m). 

Using this in ()20p . together with the change of variables i = m, j = m + n2 in the first sum, and 
i = m + rii, j = m in the second, the left-hand side of (j20p changes into 



E 



j>i>0 \ 2 

i = u {mod 2) 




i — li (mod 2) 







which allows us to combine the two sums into a single one, and so to rewrite ()20p as 

I 1, u = V = \ 



E 



i,j>0 
i = u {mod 2) 



min{i,j}\ -niin{i,j}\ 



0, otherwise, 



for all u,v>0. (22) 



Denote the above equation by E{u,v), and let c{u,v,i,j) be the coefficient of f{i,j) in E{u,v): 

/— min{i,j}\ / — min{i,j} 



c{u,v,i,j) 



u—i 
2 



.^-^•-v. 



0, 



i = u (mod 2), 
otherwise. 



(23) 



Proposition 2 Let u,v,i,j > 0. Then 

(i) c{u,v,u,v) = 1, 

(ii) c{u,v,i,j) =Oifi>uorj>v. 

Proof: Assertion (i) is clear from (j23p . To prove (ii), assume that i > u oi j > v. If i > u then 
(u — i)/2 < 0. Otherwise i < u. Then, by assumption, j > v, and so v — j — {u — i)/2 < 0. In either 
case, c{u,v,i,j) =0. □ 

Proposition [2] imphes that we can compute f{u,v) from E{u,v), provided that we have already 
computed f{i,j) for all (z,j) 7^ (n, u) such that i < n and j < v. In other words, the system (f22l) is 
a linear recurrence from which all the f{i,j) can be computed one by one, in any order compatible 
with the standard componentwise partial order on N x N. Nevertheless, we'll continue to regard 
(j22p as an infinite system of linear equations, and will rewrite it in the form Ax = b where A is 
an infinite matrix and x, b are infinite vectors. Then we can rephrase Proposition [2] as follows: 

Corollary 2 Let p : N x N ^ N be a monotonic bijection in the sense that p{a, b) < p{c, d) 
whenever a < c A b < d. Define A = [a{n, k)]^f^^Q, x = [x(A;)]^q and b = [b{n)]^^Q by 

a{n,k) = c{u,v,i,j), (24) 

x{k) = f{i,j), 

b(n) = I ^^ ^ = ^ = 1' 
1 0, otherwise, 

where c resp. f is given by [23\) resp. f21\}. {u,v) = p^^{n), and {i,j) = p^^{k). Then 
(i) A is lower-triangular with unit diagonal, 
(ii) Ax = b. 

Proof: 

(i) By Proposition [2^i), a{n,n) = c{u,v,u,v) = 1, proving that A has unit diagonal. Now 
assume that n < k. Then p{u,v) < p{i,j), hence by monotonicity of p, i > u 01 j > v. By 
Proposition [2l[|ii) , a{n,k) = c{u,v,i,j) = 0, proving that A is lower-triangular. 

(ii) Let n e N be arbitrary, and (u, f) = p^^{n). Then by (i), the sum J2T=o'^i^^^) ^i^) exists, 
and by (f23]) and (l22]) . 



J2a{n,k)x{k)=Y.c{u,v,i,j)f{z,j) = l J^ ^^^^J^J' = bin), 

k=0 i,j>0 (. 

proving that Ax = b. □ 

We can compute a particular component x{k) of the solution vector x from the finite lower- 
triangular system with unit diagonal 



where 



(k) 



W,i)]i 



i,j=0 ' 



X 






By Cramer's rule, 



x{k) 



det A 



(fc) 



det A 



(fc) 



Ak) 



det A(^) 
where A is obtained from A^ ' by replacing its last column with h^ ' . 

li k < p(l, 1) then h^^^ = and x{k) =0. li k > p(l, 1) then the last column of A^'^'' has a 
single non-zero entry, 1, at position /9(1, 1). Developing det A w.r.t. this column yields 

x{k) = detA^''^ = detT('=)detl/(*^) = det/i"^^) 

where T*- Ms a p(l, 1) x p{l, 1) lower-triangular matrix with unit diagonal, and H^ ' is the {k — 
p{l, 1)) X {k — p(l, 1)) lower-Hessenberg matrix with unit superdiagonal, composed of the elements 
in rows p{l, 1) to k and columns />(!, 1) — 1 to /c — 1 of A*- ^ Thus for Gessel numbers we have 

F(2n;0,0) = F{2n + 1;0,0) = f{2n+l,2n + l) = x{p{2n+l,2n + l)) = detii"(''(2"+^'2"+^)) 

when n > 0. For example, if we use diagonal ordering to pack the unknown f{i,j) into the vector 
X, then 

r -^ (i+J + A ^ • 
p{hJ) =1 2 J +-^' 

so p(l, 1) = 4 and /9(3, 3) = 24. Hence, if n = 1, F(2, 0, 0) = det H^'^^'^ = 2 where 
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7 Gessel numbers as multiple sums 

For any (finite or infinite) lower-triangular matrix A = [a{k,m)]k^m>o with unit diagonal entries, 
the lower-triangle elements of its inverse A~ = [a{k,m)]k,m>o are given by 

k—m j 

-a{k,m) = 5] (-IF Y. na(Ai,A,_i) (25) 

j=l m=\o<\\<---<\j=k i=l 

when k > m. Therefore for any ordering p as described in Corollary [21 

F(2n;0,0) = a(p(2n + l,2n + !),/)(!, 1)) 
where a, a, and c are given in (p5]) . (p^ . and ([23l) . respectively. 

8 The solution vector 

Here we describe a few properties of the solution vector x'"' . 

We'll think of the n^ x 1 vector x^") as consisting of the concatenation of [n^/(2n + 1)J vectors 
Uj [i = 1,2, . . .), each of length 2n+l, plus one more, of length n? mod 2n+l . Each of these vectors 
Uj, except possibly the last, consists first of a certain universal sequence of length 2i, followed by 
(2n + 1 — 2i) O's. The last one consists of as much of the next universal sequence as there is room 
for. As n increases these vectors remain unchanged, and a new one appears at the end. Each of 
these universal sequences ends in a Catalan number. 

The first several of these universal sequences are 

1, 1 

1, 2, 3, 1 

1, 5, 11, 19, 10, 2 

1, 9, 37, 85, 158, 103, 35, 5 

1, 14, 87, 332, 782, 1521, 1126, 499, 126, 14 

1, 20, 172, 911, 3343, 8004, 16056, 12941, 6765, 2296, 462, 42 

1, 27, 305, 2096, 10147, 36350, 88044, 180621, 154750, 

90681, 37178, 10254, 1716, 132 
1, 35, 501, 4300, 25927, 118472, 417565, 1020162, 

2128824, 1910006, 1217523, 570409, 193137, 44913, 6435, 429 

The sequence of next-to-last members of the above is also holonomic, but the third-from-last se- 
quence might or might not be. 

Can we find these sequences by dealing only with the corresponding sections of the matrix? 
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